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We have performed an ab initio investigation of transition metals (TMs = Mo, Ru, Co, and Pt) 
embedded in twisted bilayer graphene (tBG) layers. Our total energy results reveal that, triggered 
by the misalignment between the graphene layers. Mo and Ru atoms may form a quasi-periodic (tri¬ 
angular) array of intercalated atoms. In contrast, the formation of those structures is not expected 
for the other TMs, Co and Pt atoms. The net magnetic moment (m) of Mo and Ru atoms may be 
quenched upon intercalation, depending on the stacking region (AA or AB). For instance, we find a 
magnetic moment of 0.3 //b (1.8 /xb) for Ru atoms intercalated between the A A (AB) regions of the 
stacked twisted layers. Through simulated scanning tunneling microscopy (STM) images, we verify 
that the presence of intercalated TMs can be identified by the formation of bright (hexagonal) spots 
lying on the graphene surface. 


I. INTRODUCTION 

The synthesis of graphene and the experimental mea¬ 
surements of its outstanding properties ^ opened a path¬ 
way to explore and synthesize other single layer or few 
layers two dimensional (2D) materials such as hexag¬ 
onal boron nitride (hBN), molybdenum disulphide 
and phosphorene, a single layer of black phosphorus^. 
Furthermore the concept of van der Waals (vdW) het- 
erostructure has emerged^, layers of 2D materials, which 
interact by means of vdW forces, are piled up, aiming to 
provide a suitable set of electronic and structural prop¬ 
erties. Recently heterostructures conformed by vertically 
stacked layers of graphene/hBN/graphene have been suc¬ 
cessfully synthesized^, in these cases the different lat¬ 
tice constants and/or a relative rotation angle (RRA) 
between layer lead to moire patterns that are observed 
in scanning tunneling microscopy (STM) images. 

The synthesis of few layer graphene by chemical 
methods^’^, for instance, shows quite often a rotation 
between graphene layers with respect to the Bernal or 
AB stacking. This gives rise to a moire pattern with 
a periodicity that can be associated with the RRA be¬ 
tween layers. This structure is known as twisted bilayer 
graphene (tBG) and it has been extensively studied due 
to its peculiar electronic properties^'^^. 

In tBG, there is a so-called A A region where the stack¬ 
ing is similar to A A bilayer graphene, with one atom ex¬ 
actly above another atom of the other layer, and another 
region with AB stacking, the AB region. Between the 
AA and AB regions the stacking looks like two random 
slipped layers. The size of the tBG unit cell (UG) and of 
the AA region is increased for low RRA and the system 
can be thought as a large AA region surrounded by AB 


areas, as depicted in Fig. 1(a). Low energy electrons are 
localized in the AA region as they can not easily pene¬ 
trate the barrier formed by the AB stacking^^ The 
local density of states is significantly larger at the AA 
site for very low angles and this is precisely the region 
one observe in STM images^^’^^’^^. The electronic spec¬ 
trum also depends on the RRA. For large RRA (^ 20°), 
the system behaves as if the two layers were uncou¬ 
pled with a linear energy dispersion like in monolayer 
graphene. Below 20° there is a constant reduction of 
the Fermi velocity of charge carriers and close to 1° the 
spectrum shows zero-velocity flat bands. Scanning tun¬ 
neling spectroscopy (STS) measurements have shown an¬ 
gle dependent Van Hove singularities(vHs) near the fermi 
level^^ consistent with the behavior predicted by theoret¬ 
ical methods^’^^’^^ 

There have been recently several proposals to exploit 
the properties of graphene by intercalation of metals. An 
enhancement of the out of plane magnetic anisotropy 
have been observed in intercalated Gobalt atoms between 
graphene and Ir(lll) substratethe anisotropy also 
depends on the specific site of the unit cell where the Go 
atoms are deposited In some cases the intercalated 
Go atoms forms a quasi-periodic array of trimers between 
graphene adlayer and the SiG(OOOl) surfaceBesides, 
the geometry of adatoms adsorbed on the graphene sur¬ 
face, might play a very important role on the formation of 
quantum spin and anomalous Hall phases in graphene . 

In this work we propose to take advantage of the moire 
structure in twisted bilayer graphene to create a peri¬ 
odic array of intercalated transition metal(TM) atoms. 
If there is any energetic preference for the atoms to be lo¬ 
calized in a given region between the layers of the twisted 
unit cell it might be possible to control the distance be- 
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tween them by manipulating the RRA between layers, 
allowing or not the interaction between atoms and con¬ 
trolling the electronic/magnetic properties of the system. 

We examine, based on ab initio total energy calcula¬ 
tions, the energetic stability and the electronic properties 
of TMs (TM=Mo, Ru, Co, and Pt) intercalated between 
tBG layers. We found that Mo and Ru atoms exhibit an 
energetic preference for the A A stacking region (TMaa) 
of the twisted bilayer unit cell, while no such a preference 
is found for Co and Pt. This suggest that Mo and Ru 
intercalated atoms might form a quasi-periodic triangu¬ 
lar array embedded in the twisted layers. At the same 
time the net magnetic moment of Mo, Ru and Co atoms 
are somewhat quenched upon their intercalation, for in¬ 
stance, we found a magnetic moment of 0.3 /xb foi* Ru 
adatoms intercalated in the AA stacking region. Through 
STM simulation, we show that there is an increase of the 
electronic density of states on the (graphene) C atoms 
bonded to the TM incorporation sites. 


II. METHOD 

A. Geometry of the Unit cell. 

The system we studied is composed of two coupled 
graphene layers with a given RRA between them. We 
built a commensurate unit cell following a procedure ex¬ 
posed in several previous works^^’^^. Starting from a 
stacked AB bilayer graphene one rotate a lattice point to 
an equivalent location, a vector r = mai + na 2 is rotated 
to t = ndi + 17102, where di and 02 are the graphene 
bilayer lattice vectors; n^m are integers. These integers 
are the only information needed to completely character¬ 
ize the unit cell, the rotating angle, the unit cell vectors, 
and the numbers of atoms in the UC can be expressed in 
terms of them^^’^^. The UC constructed this way con¬ 
tains one site, known as A A, where one atom is exactly 
above another atom of the other layer, there is another 
site, AB, where an atom of the top layer is exactly in the 
center of a lower layer hexagon, we have also a BA site, 
equivalent to the AB site but with the difference that an 
atom in the lower layer is in the center of an hexagon 
of the upper layer, there is also another area of the UC 
where the stacking is none of the above mentioned we 
call it slipped AB (SAB), see Fig. 1(a). As the UC gets 
larger, very small RRA, the number of atoms with almost 
exact AA or AB character increases^^ creating large re¬ 
gions with almost the same stacking. The distance be¬ 
tween two AA regions (periodicity of the superlattice) 
can be calculated by the equation D = ao/2 x sin(^/2), 
where oq = 2.46 A and 0 is the RRA between layers. 


B. DFT calculations 

The calculations were performed within the density 
functional theory (DFT) approach, as implemented in 
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FIG. 1. (a) Structural model of twisted graphene bilayer with 
RRA of 9.34° (tGB-148). AA, AB, BA and SAB stacking 
regions are indicated (see text) and the arrows represent the 
supercell lattice vectors, (b) tBG-148 total electronic density 
of states (DOS) and (c) projected DOS (PDOS) on the carbon 
atoms near the AB (Cab) and AA (Caa) stacking regions. 
In (b) and (c), we have considered a gaussian broadening of 
0.1 eV in the discrete electronic energies, and the Fermi levels 
were set to zero. 


the Quantum-Espresso package^^. The Kohn-Sham or¬ 
bitals were expanded in a plane-wave basis set, with an 
energy cutoff of 35 Ry. The twisted graphene bilayers 
were simulated considering supercells with 148 (tGB- 
148) and 244 (tGB-244) atoms, corresponding to RRA 
{0) of 9.34° and 7.34°, respectively, and vacuum regions 
of 14 A. calculations. The atomic positions were fully 
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relaxed, by including van der Waals (vdW) interactions 
within the self-consistent vdW-DF approach as described 
in Refs.^^~^^. For the structural relaxations we have con¬ 
sidered a force convergence tolerance of 15 meV/A and 
the Brillouin zone (BZ) was sampled by 4x4x1 (2x2x1) 
k-points in the case of tGB-148 (tGB-244). After the 
structural relaxations, spin polarized calculations were 
performed with 8x8x1 (5x5x1) k-points BZ sam¬ 
pling in the case of tGB-148 (tGB-244). 


III. RESULTS 

We begin examining the electronic and structural prop¬ 
erties of the twisted bilayer graphene crystal in compar¬ 
ison with recent experimental and theoretical findings. 
In Fig. 1(a) we show the structural model of tGB-148, 
indicating the AA and AB stacking regions. At equilib¬ 
rium, the averaged interlayer distance (Az) is 3.679 A; 
slightly larger than the distance obtained for perfect bi¬ 
layer graphene (AB or Bernal stacking). There are no 
chemical bonds between the graphene sheets, the inter¬ 
layer interaction is governed by vdW forces, we found a 
binding energy of 27meV/G-atom very close to the ex¬ 
perimental value of graphite^^’^^. In the case of tGB-148, 
RRA = 9.34°, a slight renormalization of the Fermi ve¬ 
locity and the presence of vHs, associated with the M 
point of the BZ are expected^^. The energy of these 
vHs depends on the RRA, diminishing as the angle gets 
lower found that the vHs lie at ±0.64 eV with 

respect to the Dirac point, as shown in Fig. 1(b) in agree¬ 
ment with previous experimental results^^. 

The moire structure gives rise to a periodic potential 
that have a huge impact on the electronic dispersion and 
at the same time promotes a charge density redistribu¬ 
tion confining low energy electrons to the AA region, pro¬ 
ducing the bright spot observed in STM experiments on 
graphite or tBG 

The origin of these bright spots were the subject of 
a debate not so long ago^^. Two theories were vindi¬ 
cated: (i) there is an increase of the electronic density 
of states around the AA region, when compared with the 
AB area and (ii) the AA and AB regions present different 
equilibrium geometries, and the interlayer distances Az 
(topographic corrugation) are different. Both arguments 
are valid and indeed, (i) can be verified in Fig. 1(c), by 
projecting the electronic DOS on the G atoms near the 
A A (Gaa) and AB (Gab) stacking regions. The elec¬ 
tronic density of states is larger on Gaa- la addition, 
comparing the equilibrium geometries, we find that the 
interlayer distance in the AA region (Azaa = 3.692 A) is 
larger by 0.04 A than in the AB one (Azab = 3.653 A), 
as proposed in (ii). By reducing the twist-angle to 7.34° 
(tGB-244), the lateral distance between the AA and AB 
regions increases, when compared with tGB-148, and we 
find practically the same (averaged) interlayer distance, 
3.681 A. However, the interlayer distance in A A region 
increases to Azaa = 3.765 A, while in the AB region 


we find practically the same distance (Azab = 3.650 A). 
These small changes in the interlayer distances suggest 
that the origin of the STM images is mainly electronic 
rather than geometrical^^. 

We found pertinent to investigate whether this behav¬ 
ior facilitates the creation of periodic arrays of interca¬ 
lated atoms between twisted graphene layers. To such 
end we studied the energetic stability of transition metal 
atoms (TMs: Mo, Ru, Go, and Pt) intercalated at differ¬ 
ent sites of the tBG unit cell. A measure of the energetic 
stability is the binding energy of the TMs embedded in 
the twisted systems, the binding energy (E^) can be ob¬ 
tained from: 

= E[tBG] ± E[TM] - F;[TM/tBG]. 

Where F^[TM/tBG] is the total energy of the final system, 
a tBG intercalated by a TM atom. F^[tBG] and F^[TM] 
are the total energies of the separated components, a tBG 
and an isolated TM atom, respectively. 

We calculated the binding energy at several locations 
within the UG and found that for Go and Pt the differ¬ 
ences are small, there is not any clear preference for any 
region. However for Mo and Ru the difference can be 
as large as 1 eV, and the highest binding energy is ob¬ 
tained in the AA region and the lowest in the AB one. In 
the rest of the UG the results are in between these two 
extreme values. 

In addition we have calculated the binding energy but 
with the TM atoms adsorbed on the surface site (an atom 
on top of the upper layer). We obtained, for the tBG-148 
structure, E^ of 3.48, and 2.47 for Mo and Ru adatoms 
lying on the hollow site of the AA region, and 1.44 eV for 
Pt adatom on the bridge site. These values are signifi¬ 
cantly smaller than the ones obtained for the intercalated 
atoms. 


TABLE I. Binding energies (E^ in eV/TM-atom) of TMs 
intercalated in the AA and AB stacking regions of twisted 
graphene bilayer, tBG-148 and tBG-244. 


TM 

e = 9.34° 
AA 

(tBG-148) 

AB 

0 = 7.34° 
AA 

(tBG-244) 

AB 

Mo 

5.79 

4.64 

5.88 

4.73 

Ru 

3.76 

3.16 

3.80 

3.20 

Co 

2.74 

2.73 

2.74 

2.73 

Pt 

1.66 

1.63 

1.77 

1.71 


Table I shows the values of the binding energy for the 
AA and AB regions for all the intercalated TM con¬ 
sidered. There is an energetic preference (higher bind¬ 
ing energies) for Ru and Mo atoms intercalated in the 
AA region. For the tBG-148 system, the Moaa(Ruaa) 
structure is more stable than Moab(Puab) by 1.15 eV 
(0.60 eV). Such energetic preference has been maintained 
in tBG-244. That is, the binding energies of Moaa 
and Ruaa are higher by 1.15 and 0.60 eV, respectively, 
compared with the ones of Moab and Ruab- On the 
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other hand, the energetic preference for the A A stack¬ 
ing region reduces to 0.03 eV (0.06 eV) for PtAA/tBG- 
148 (PtAA/tBG-244), while the Goaa and Goab systems 
present practically the same binding energy. In the ap¬ 
pendix we show the results of the binding energies for Ru 
and Mo in other regions of the UG. In general within the 
AA region the differences are small (around 1%) and in 
the SAB and AB regions the binding energies are around 
5% and 20 % lower, respectively, than in the AA region. 

In Figs. 2(a) and 2(b) we show the structural mod¬ 
els of Mo atoms intercalated in the regions A A and AB, 
the structural models of Ru and Go are very similar. In 
the AA stacking region the TM atom lies on the center 
of one of the empty hexagons (Fig. 2(a)), in this case 
the TM atom is bonded to 12 carbon atoms^^. On the 
other hand, in the AB stacking region (Fig. 2(b)), the 
TM atom lies on top of one atom of one layer, and on 
the hollow site of the other one. In this case, the co¬ 
ordination of the TM atom is reduced to seven. The 
interlayer distance, upon the presence of intercalated Mo 
and Ru atoms, Az = 3.73 A, is slightly larger compared 
with the one obtained for the pristine twisted system. 
Whereas, we find Az around 3.85 A for Moab and Ruab, 
thus suggesting that those structures induce a larger local 
strain when compared with their Moaa and Ruaa coun¬ 
terparts. Differently from the other TM atoms, the Pt 
atoms lie on the bridge site of the graphene layer in both 
stacking regions. Figs.2(c) and 2(d). Thus in both cases 
Pt atoms are bonded to four carbon atoms. Mo and Ru 
atoms, as d-orbital metals, tend to have as many bonds 
as possible when they are in contact with organic mate¬ 
rials while Pt atoms behave more like noble elements. 

For Go atoms the different coordination in the AA and 
AB stacking region do not lead to large differences in 
the total energies. This is due to the fact that Go co¬ 
valent radius is smaller than Mo and Ru radius, thus it 
is expected because the covalent bonds between Go and 
carbon atoms are not as strong as Mo(Ru)-G bonds. 

We proceed now to quantify the local strain in the 
twisted graphene bilayer, upon the intercalation of the 
TM atoms, we calculate the deformation energy (F’^®-^), 

^def _ ^strain _ ^prist 

here and represent the total energies of 

pristine (full relaxed) and strained tBG-148 system, re¬ 
spectively. The latter term was obtained by considering 
the atomic positions of the deformed (due to the presence 
of the intercalated TM) tBG-148 structure. The defor¬ 
mation energies for Mo and Ru atoms intercalated in the 
A A stacking region Moaa and Ruaa are 0.12 and 0.09 eV 
respectively, while it increases by about O.IeV = 

0.24 (0.19) eV) for Moab (Bclab)- That is, the local 
strain also contributes to the energetic preference for the 
AA stacking regions for these two transition metals. On 
the other hand, for intercalated Co and Pt atoms, the 
differences are smaller, 0.024 and 0.015 eV, respectively. 

Summarizing the total energy results; we have found 
that the most stable place for Mo and Ru atoms, in¬ 


tercalated in tBG, is in the A A stacking region. This 
total energy difference may be explained by the differ¬ 
ent chemical coordination and also by the strain energy 
differences. 

Therefore our results indicate that Mo and Ru atoms 
may form a quasi-periodic triangular array when they are 
intercalated in twisted graphene layers, the most stable 
position is within the AA region and the energy differ¬ 
ences with other regions are larger than room tempera¬ 
ture thermal energy. The periodicity of these arrays can 
be tuned by the twist angle. 

The relative charge concentration around the AA re¬ 
gion increases as the angle gets lower we expect then that 
the binding energy at the AA region will be larger when 
we compare with the rest of the UC for lower angles as 
the charge density should be proportional to the strength 
of the covalent bond between TM and Carbon. 



FIG. 2. Equilibrium geometry of MoAv/tBG, top-view (al) 
and side-view (a2), and Moab /tBG TM(AB), top-view (bl) 
and side-view (b2). The same structural model has been ob¬ 
tained for Ru/tBG and Co/tBG. Equilibrium geometry of 
PtAA/tBG (c) and PtAB/tBG (d). 

We will focus now on the electronic and magnetic prop¬ 
erties of intercalated TM atoms. Mo and Ru have differ¬ 
ent magnetic moments (m) when we compare with the 
adsorbed ones, as can be noticed in Table II. The mag¬ 
netic moment of Mo is quenched when it is intercalated 
in both A A and AB stacking regions. Ru magnetic mo¬ 
ment is also quenched in the more energetically stable 
AA stacking region but not in the AB region. Those 
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FIG. 3. (left) Total density of states of states (DOS), shaded 
plots, and the projected density of states (PDOS), solid lines, 
of the TM-4d (z^, xz, yz, xy, and x^ — y^) orbital, (right) pro¬ 
jected density of states of C-2pz orbitals of C atoms of the 
graphene layer A (see Fig. 2) bonded to the TM (red solid 
line) and far from the TM (black solid line); for Mo/tBG-148 
(a)-(d), and Ru/tBG-148 (e)-(h). The (spin polarized) DOS 
and PDOS are in arbitrary units, and the Fermi level was set 
to zero. 


changes on the magnetic moment are ruled by the differ¬ 
ent TM hybridization with the tBG host. The magnetic 
moments of Co is almost unchanged when comparing ad¬ 
sorbed and intercalated results. This is another evidence 
for the fact, mentioned above, that Co-carbon bonds are 
weaker in comparison with Mo-carbon and Ru-carbon 
bonds. 

There is an electronic charge transfer from the Mo 
and Ru atoms in the A A region to the graphene layers. 
This lead to n-type doped tBG, with charge transfers 


TABLE II. Magnetic moment (m in /xb) of the TMs adsorbed 
(mads), and intercalated in the AA and AB sites (mAA and 
hiab) of the twisted (tBG-148) graphene layers. 


TM 

rUads 

mAA 

mAB 

Mo 

2.0 

0.1 

0.0 

Ru 

1.9 

0.3 

1.8 

Go 

1.2 

1.4 

1.2 

Pt 

0.0 

0.0 

0.0 


of 3.2 and 5.1 X lO^^cm for Mo and Ru atoms, re¬ 
spectively. Also, the local hybridization of TMs in AA 
and AB regions are different, as shown by the calculated 
spin-polarized electronic density of states (DOS), shaded 
plots in Fig. 3, and projected DOS (PDOS, solid lines in 
Fig. 3). The Mo-4d electronic states are mostly localized 
below the Fermi level {Ep — leV), and the net mag¬ 
netic moment of 0.1 /ip comes from a small occupation 
unbalance of Mo-4d22 orbitals, indicated by green solid 
line in Fig. 3(a). It is worth noting that although the 
symmetry of the intercalated Moaa atom is slightly 
missed, due to the twist angle, the double degeneracy of 
4dxy/4dx2_y2 and 4dxz/4dyz (levels e and e') have been 
preserved. In addition, as presented in Fig. 3(b), there is 
an increase of the DOS projected on the C atoms bonded 
to Moaa, Caa in Fig- 2(a2). 

The hybridization of Mo atoms in the AB region with 
the graphene layers has been strengthened due to the for¬ 
mation of Mo-C bond with a C atom of the graphene layer 
B [Fig. 2(b2)]. In this case, the occupation of Mo-4dz2 
reduces by around 20%, followed by an increase of the 
PDOS on the C atoms bonded to Moab (Cab), Figs. 3(c) 
and 3(d). 

In Figs. 3(e) and 3(f) we present the DOS and PDOS 
of RuAA/tBG-148. Here, the double degeneracy of the 
levels e and e', composed by Ru-4dxy/4dx2_y2 and Ru- 
4cixz/4dyz have been maintained [Fig. 3(e)], while the Ru- 
4dz2 (lying at Ep —1.9 eV) weakly interact with the tBG- 
148 host. The net magnetic moment of 0.3 /xp comes 
from the partial occupation unbalance of Ru-4dxz and 
4dyz orbitals, as well as from the 2pz orbitals of Caa 
atoms. Fig. 3(f). Similarly to the MoAp/tBG, the lo¬ 
calized character of RuAA-4dz2 orbitals has been missed 
in RuAp/tBG, Fig. 3(g). In addition, we find that spin- 
down states, composed by Ru-4dxz and -dd^^ orbitals, be¬ 
comes unoccupied increasing the net magnetic moment to 
1.8/xp. The electronic DOS of the carbon atoms bonded 
to Ruap, Cap, increases compared with the other C 
atoms of the graphene layer [Fig. 3(h)]; contributing to 
the net magnetic moment of Ruap/IBG. 

The changes in the electronic density of states on the 
graphene sheet, induced by the intercalated TMs, can be 
visualized through STM images. Here, based on Tersoff- 
Hamman approach, we simulate a STM image of occu¬ 
pied states, within Ep — 1 eV, of Moaa/ and RuAA/tBG- 
148. As depicted in Figs. 4(a) and 4(b), we find bright 
hexagonal spots lying on the C atoms above the inter- 
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calated Moaa and Ruaa atoms, respectively. Those 
spots form a triangular array, and the brightness can be 
assigned to the larger PDOS of Caa atoms [Figs. 3(b) 
and 3(f)]. In addition to the electronic contribution, 
the graphene surface corrugation may play an impor¬ 
tant role. As discussed above, the interlayer distance 
increases slightly upon the presence of intercalated TMs, 
e.g. Az = 3.69 ^ 3.73 A for Mo and Ru in the AA 
stacking region. 


(a) 

(b) 

(c) 

(d) 


FIG. 4. Simulated STM images of occupied states {Ep — 1 eV) 
of a tBG-148 system with iutercalated Moaa (a) aud Ruaa 
atoms (b). Simulated STM images of the same system with¬ 
out the iutercalated atoms but with the corrugatiou iuduced 
by them, Moaa (c) aud Ruaa (d). 


In order to verify the importance of such a corruga¬ 
tion on the STM images, we simulated STM images of 
the twisted graphene (tBG-148) by removing the inter¬ 
calated atoms; but keeping the equilibrium geometry of 
the Moaa/ and RuAA/tBG-148 systems. Our simulated 
STM results presented in Figs. 4(c) and 4(d) show that 
indeed the bright spots follow the topographic corruga¬ 
tion of the graphene surface. 

Several studies have been addressed to the forma¬ 
tion of magnetic clusters at the graphene/surface inter¬ 
face In those systems, the formation of (self- 

assembled) periodic structures is ruled by a suitable syn¬ 
ergy between the hexagonal lattice of the graphene layer, 
and the atomic geometry of the surface. Our results sug¬ 
gest that the array of Mo and Ru atoms intercalated in 
the twisted graphene layers can act as a seed to the for¬ 
mation of these type of clusters. Indeed, our preliminary 
results indicate that the formation of Ru-dimers inter¬ 
calated in the AA stacking region is quite likely. For 
Ru-dimers we find of 4.00eV/atom, whereas for Ru- 
trimers we obtained 4.05eV/atom. Further studies re¬ 
garding the formation of TM clusters in twisted graphene 
layers are in progress. 


IV. SUMMARY 

We have performed an ab initio investigation of TM 
atoms (Mo, Ru, Co and Pt) intercalated in twisted 
graphene layers. Our results show that Mo and Ru are 
significantly more stable energetically at tBG A A stack¬ 
ing region in comparison with the rest of the unit cell. 
For Co and Pt atoms the energy differences between re¬ 
gions are not relevant. These results indicate that it is 
possible a formation of a triangular array of Mo and Ru 
atoms, ruled by the twist angle between the graphene 
layers. In contrast, formation of such a structure is not 
expected for Co and Pt atoms. The energetic stability 
of Mo and Ru atoms at tBG AA stacking regions can be 
explained in terms of the coordination number of the TM 
atoms with carbon atoms and also differences in strain 
energy. The net magnetic moments of Mo is quenched 
upon its intercalation in A A or AB regions, while Ru 
magnetic moment is 1.9 in the adsorbed system and goes 
to 0.3 /iB when intercalated in the A A region. There is 
an increase of the electronic density of states, near the 
Fermi level, on the graphene carbon atoms bonded to the 
TM atoms, identified through STM images. We find the 
formation of a triangular array of hexagonal bright spots 
lying on the C atoms nearest neighbour to the interca¬ 
lated Mo and Ru atoms. These spots in the STM images 
may be used to experimentally identify the intercalated 
atoms positions. 

The authors acknowledge financial support from the 
Brazilian agencies CNPq, CAPES and FAPEMIG, the 
computational support from CENAPAD/SP and Chilean 
EONDECYT grant 1130950. 

V. APPENDIX 

In Eig. 5 we show other positions of the intercalated 
TM atoms inside the tBG unit cell(tGB-148). The results 
are for Ru and Mo intercalated atoms since these are 
the systems with large energy differences between AA 
and AB regions. In these cases in order to preserve the 
RRA unchanged during the atomic relaxation process, 
the C atoms of graphene layers are allowed to relax only 
along the z-direction; such a constrain is not applied to 
the intercalated TMs. In Table III we present the total 
energy differences (AE) with respect to the energetically 
more stable TMaa position, 

ae; = e;[tMaa]-^[tMx], 

where £’[TMaa] represents the total energy when the 
TM is located at the A A region, and E[TMx] is the total 
energy of the other locations (X = AA-i, SAB, SAB-i, 
AB, and AB) depicted in Eig. 5. These results indicate 
that indeed, there is an energetic preference for the AA 
sites. The total energy difference (AE) increases from 
TMaa to TMab- The results of Tables I and III lead 
us to conclude that the in plane relaxation {x and y di¬ 
rection) plays a minor role on the energetic preference of 
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TM atoms. In Table I, where in plane relaxation is con¬ 
sidered, the results are very close to the ones obtained in 
Table III where it has not been allowed. 



FIG. 5. Structural models of TM/tGB-148, for TM=Ru and 
Mo between the energetically more (less) stable TMaa/TGB- 
148 (TMab/TGB-148) configurations. 


TABLE III. . Total energy difference, AE = E[TMx] — 
E[TMaa] (in eV) between the energetically more stable 
TMaa/TGB-148 (TM = Ru and Mo), and the TMx ones 
presented in Fig. 5. 


TMx 

Ru 

Mo 

TMaa 

0.000 

0.000 

TMaa-i 

0.085 

0.120 

TMsab 

0.168 

0.249 

TMsAB-i 

0.290 

0.541 

TMab-i 

0.296 

0.478 

TMab 

0.583 

1.185 



FIG. 6. Structural models of TM/tGB-148, for TM=Ru and 
Mo atoms embedded in AA-like regions, AA and AA'. 

In addition, we compare the total energies for Rn/ and 
Mo/tGB-148 embedded in AA-like regions, indicated as 
A A and A A' in Fig. 6. The atomic positions of the TM 
and tGB are free to relax. Those AA and AA' configura¬ 
tion are very close in energy; we find differences of 0.007 
and 0.024 eV. 
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